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Possibility of structured mixed phases at first order phase transitions is exam- 
ined with taking into account of charge screening and surface effects. Hadron- 

' quark phase transition in dense neutron star interior is considered, as concrete 

O ' 

example. 
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Ref. [1] suggested presence of a wide region of mixed phase at any first order phase 
transitions in mult i- component systems of charged particles. Existence of structured mixed 
phase in dense neutron star interiors would have important consequences for equation of 
state, also affecting neutrino emissivities [2], glitch phenomena and r modes, cf. [3,4]. How- 
ever inhomogeneity effects of the field profiles were disregarded in these treatments. On the 

X" 

other hand, ref. [5] demonstrated that for the appearance of the structured mixed phase 
the Coulomb plus surface energy per droplet of the new phase should have a minimum, as 
function of the droplet radius. Corrections to the Coulomb solutions due to screening effects 
were disregarded. In spite of the question has been rised long ago, cf. [1-6], up to now there 
is no consistent treatment of these effects. Therefore, further study of the screening and 
surface effects seems to be of prime importance. 

Consider the structured mixed phase consisting of two phases I and II. We suppose the 
lattice of droplets (phase I) each placed in the Wigner-Seitz cell (the exterior of droplets is 
phase II). Each droplet in the cell occupies the domain D l of volume v l separated by a sharp 
boundary 3D from matter in phase II ( region D u of volume v 11 ) . We exploit thermodynamic 
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potential (effective energy) per cell composed of a density functional [7], 

n = E\p\- /4 J di dfp\ - J du drpY, (i) 

E[p] is the energy of the cell, p = {p\, p\ 1 } are densities of different particle species, i = 
1, iV 1 in phase I and % — 1, N n in phase II, N 1 , N u are total number of particle species 
per cell in phases I and II. Summation over the repeated Latin indices is implied. Chemical 
potentials p], pf are constants, if, as we assume, each phase is in the ground state. We also 
assume that matter in phase I or II is in chemical equilibrium by means of the weak and 
strong interactions. Equations of motion, |pr = 0, render 

tf=^d, a = {I,n>. (2) 

The energy of the cell consists of four contributions: 

E[p] = drei^M + J m dfe^ n+str [pf} + dSe s [p] + E v . (3) 

The first two contributions are the sums of the kinetic and strong-interaction energies and 
es[p\ is the surface energy density, which depends on all the particle densities at the boundary 
dD. One may approximate it, as we shall do, in terms of surface tension a. 1 Ey is the 
Coulomb interaction energy, 

E V = 1 - l drd? QiPi ™i P f\ (4) 
2 J | r — r' | 

with Qi being the particle charge (Q = — e < for the electron). 
Then equations of motion (2) can be re-written as 

f) f Ot 

tf = ^f L -N?> a (V a -V°), N?> a = Q?/e, (5) 



1 In reality there is no sharp boundary between phases and all the densities and constant chemical 
potentials are defined in the whole space (c.f. models of atomic nuclei). Then particle densities, 
pf, are changed continuously. The surface energy is given by integration of e a over a narrow region 
around surface, where pf change sharply. Our formulation is also applicable to this case. 



where V a (r) is the electric potential generated by the particle distributions, which can be 
shifted by an arbitrary constant (V°) due to the gauge transformation, V — > V — V°, 

v > ^-/^£^) + v-«J m ?eB ' (6) 

J \r-r'\ [V"(r), feD" 

Formally varying eq. (5) with respect to V a or p" we have the matrix form relation, 

do a 

A a Ll — /v ch ' a A a R a — ft-, (7) 
where matrices A and B are defined as 

\2 jpa 



^ -^kin+str r>a _ ^_P_ 

Eqs. (7), (8) reproduce gauge-invariance relation, 



/|« _ kin+str _ /n\ 

ij ~ 5p?5p° ' ij ~ dpf { } 



dpi = (9) 

clearly showing that const ant- shift of the chemical potential is compensated by gauge trans- 
formation of V a . Hence chemical potential p" acquires physical meaning only after fixing 
of the gauge of V a . 

Applying Laplacian (A) to the l.h.s. of eq. (6) we recover the Poisson equation (r£ D a ), 

AV a = 47re 2 p ch ' a = 4vreQ> J Q . (10) 

The charge density p ch ' a as a function of V a is determined by equations of motion (5), so 
that eq. (10) is a nonlinear differential equation for V a . The boundary conditions are 

V l = V n , W r = W n , re 3D, (11) 

where we neglected a small contribution of surface charge accumulated at the interface of 
the phases. We also impose condition W u = at the boundary of the Wigner-Seitz cell, 
which implies that each cell must be charge- neutral. Once eqs. (10) are solved giving V a 
and the potentials are matched at the boundary, we have density distributions of particles 
in the domain D a . 



Note that there are two conservation laws relevant in neutron star matter: baryon num- 
ber and charge conservation. These quantities are well defined over the whole space, not 
restricted to each domain. Accordingly the baryon number and charge chemical potentials 
(Hb and /iq) , being linear combinations of /if, become constants over the whole space, 

=^b = Vb, fJ-Q = = Hq- (12) 

This fact requires two conditions for fif at the boundary dD, which prescribe the conversion 
manner of particle species of two phases at the interface. In other words, charge and baryon 
number densities should be continuous across the boundary due to eq. (5). 2 In particular, 
electron chemical potential is equal to the charge chemical potential (/xq = //") and its 
number density is related as 

pa _ (gg -v a + yy (13) 

from eq. (5). Note that this is a gauge- invariant quantity. 

Once eq. (5) is satisfied, pressure becomes constant in each domain, 

-P°v° = df {e^ n+stI [ P ?} - l^ptV a - tfp?} . (14) 

Hence, the extremum condition for Q with respect to a modification of the boundary of 
arbitrary shape (under the total volume of the Wigner-Seitz cell being fixed) reads 

P- = P- + ^, (15) 

S is the area of the boundary dD and a is surface tension. The boundary of the cell does 
not contribute since all the densities are continuous quantities at this point. Eq. (15) is 



2 Each particle density is not necessarily continuous across the boundary, since it is only defined 
in each phase, while densities of leptons are well defined over the whole space. When particles of 
the same species i are allocated in both domains and the conversion of particle species becomes 
trivial, we must further impose the relations, fi\ = n] 1 , and p\ = p\ l at the boundary. 
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the pressure equilibrium condition between two phases. 3 Thus we satisfy Gibbs conditions 
(12), (15) in our formalism. 

The Debye screening parameter is determined by the Poisson equation if one expands 
the charge density in 5V a = V a — V" around a reference value V". Then eq. (10) renders 



A5V a = Ane 2 p ch ' a (V a = V r a ) + {n a {V a = V r a )) 2 5V a 

d 2 p ch,a~ 



(16) 



+27re" 



(dV a f 



(sv a y + ..., 



ya—ya 



[ K a (V a = V T a )Y = Ane 



dp' 



dV 



= 4 TrQtQ 



1 ^3 f),,a 



(17) 



ya—ya 



where we used eq. (9). Then we calculate contribution to the thermodynamic potential 
(effective energy) of the cell up to 0((5V a ) 2 ). A proper electric field energy of the cell is 



E V - f drV + f dre" - [ ^ldr + [ (W " )2 df 
ti v - J Di dre v + J^drey- ^ dr + J^ ^ dr; 



(18) 



that in the case of unscreened distributions is usually called the Coulomb energy. Besides 
the terms given by (18), there are another contributions arising from effects associated 
with inhomogeneity of the electric potential profile, due to implicit dependence of partial 
contributions to the particle densities on V^ u \ We will call them correlation terms. Then 
taking pi as function of V a we expand e£ in+str in 5V a : 

dpi 



e kin+str[p] — e kin+str(Pr(^ a )) + 



(pf + Nf ha (V a - V )) 



1 

+ 2 



a\2 



(««) 



+ (pf + Nf l ' a (V a -V' 



d 2 p a 



dV a 

a\2 



5V a 



ya — ya 



(sv a y + 



(19) 



ya — ya 



Ane 2 ^ 1 v '> (dV a ) 2 _ 

We used eqs. (2), (6), (7) and (8) in this derivation. The term in (19) oc (V a - V°)(5V a ) 2 is 
actually of the higher order of smallness than other terms if V a ~ Vq ~ V T a . Using expansion 



3 As noted in footnote 1, in realistic problem with continuous density distributions, i.e. in absence 
of sharp boundary, the contribution of the surface energy is absorbed into P a . Hence P l = P 11 in 
such a more detailed treatment. 
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-tip? 



-[/« 



a dp° 



dV r - 



5V a - 



a d 2 P a 



P 



(sv a y + . 



(20) 



we obtain the corresponding correlation contribution to the thermodynamic potential f2 cor = 



cor vi 



kin+str 



(p?(v r a )) - pTpUK") - p ch > a (v r a )(v r a - v°) 



^V r a -V°)AV a | (K a (V T a )f(5V a f | 



(21) 



47re 2 87re 2 
where we also used eqs. (16) and (17). In general V} ^ V^ 1 and they may depend on the 
droplet size. Their choice should provide appropriate convergence of the above expansion in 
8V. Taking = V* 1 = V T = const we find 

{^{V T )f{V a -V T f 



87re 2 



+ const, 



(22) 



and one may count the potential from the corresponding constant value. 

In the following we consider the hadron-quark phase transition, as an example. We 
suppose the lattice of spherical droplets of the radius R placed in the Wigner-Seitz cell 
of the radius i?w We also assume the quark matter inside the droplet, as phase I, and 
the hadronic matter outside in the cell, as phase II, both divided by a sharp boundary 
r = R. The quark matter consists of u, d, s quarks and electrons and the kinetic plus strong 
interaction energy density is given by 



^kin+str 



3tt 2 / 3 



1 



b/3 , 4/3 4/3 , ™IpT 
37l ) Pu +Pd +Ps + ^4/3 



+ B + 



(37T 2 Pe ) 4 / 3 

4tt 2 ; 



where B is the bag constant, a c is the QCD coupling constant, m s is the mass of strange 
quark. Last term is the kinetic energy of electrons. 

The hadronic matter consists of protons, neutrons and electrons and the kinetic plus 
strong energy density is given by 



4Utr ^ 4 m [Pn] + eflpj + e pot [p n ,p p ] + (3?r ^ e 2 )4/3 ; 



(23) 



where eY n [pi],i = n,p are standard relativistic kinetic energies of nucleons, while e pot is the 
potential energy contribution we take here in the form 
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Po -l-o \ Po / 

+ Qat(p„ + P P )(^^-l), (24) 

p is the nuclear density (p — 0.16fm~ 3 ) and constants ebmd, K , C sat are determined to 
satisfy the nuclear saturation properties. 

We use chemical equilibrium conditions for the reactions u + e <-> s, d <-> s, and n <-> p + e 
in each phase, 

Pn - Ps + Pe = 0, p d = p s , p„ = p p + p e , (25) 

and the conversion relation at the boundary, 

Ps = Pn = 2p d + p u , (26) 

which yield relations between quark and nucleon chemical potentials. 4 Using these condi- 
tions we obtain, cf. [5], 

p ch,i ^ (l - 

where the electron contribution is omitted as small (~ (V 1 ) 3 ) and we fixed the gauge by 
taking V° = -fi Q = -p e . 

Poisson equation (16) with p 011 ' 1 from (27) describing electric potential of the quark 
droplet can be solved analytically. For r < R with the boundary condition \V l {r — > 0)| < oo, 
we find 

V 1 = sh^V) + Ul (28) 

with an arbitrary constant V$. For the Debye parameter n l and for the constant U\ we 
obtain: 



4 Other conversion relation p p = 2p u + p^ is then automatically satisfied. 



2p|^ 
9tt 2 



1 + 



fi B m 2 s 




(27) 
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(K , )2 ^ (l _2 ?)i v , „ (29) 

Note that solution (28) is independent of the reference value V} in this case, c.f. (17), since 
pc\\ j n ^27) is linear function of V 1 in approximation used. 

For phase II, expanding the charge density p ch,n around a reference value, p ch ' 11 ~ 
p p (V u = V r 11 ) + 5p p - p e (V u = V v 11 ) - 5p e , and using eqs. (7), (8), (17) and (23) we find up 
to linear order 

5 Pp ^C 1 (V u (r)-V r u ), 5p e = { ^(V u (r)-V r 11 ) (30) 
Co = ^ - ,f 2 + ^, Pfp = ^(V- = 0) 1/s - 
For r > R, the Poisson equation with the boundary condition V'\ Rw = yields 

V u = V^jch (n l \r - R w )) (1-5) + U* 1 (31) 
with an arbitrary constant V U , where the constant Uq 1 is given by 

4 2 Ch,ll/T/Il _ T/II\ 

= - — + ^ < 32 > 

We will further drop numerically small 8 correction. We take the reference value V} 1 = 
V(r — Rw)- Since size of the Wigner-Seitz cell Rw is substantially larger than R, we 
have V(r = Rw) — V hnlk , where V hnlk is constant bulk solution of the Poisson equation, 
(\/ bulk ) 3 = -pi = -37T 2 p p (V rbulk = -p e ), that coincides with the local charge-neutrality 
condition for the case of the spatially homogeneous matter. Hence we find Uq 1 ~ V^ 1 ~ — p e . 
5 The charge screening in the external region is determined by the Debye parameter 

= 4eV| + 4A 

7T G 

where second term is contribution of proton screening. Taking p 1 ^ = 1.5po, p e — 170 MeV, 
Pb = Pn — 1020 MeV, a c ~ 0.4, we roughly estimate typical Debye screening lengths as 



5 Thus solutions (28) and (31) are proved to be consistent with eq. (6) for the gauge choice, 

V° = -Me- 
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A^> = 1/k 1 ~ 3.4/mTr, and A" = ~ 4.2/m^, whereas one would have A" ~ 8.5/m w , if 
the proton contribution to the screening (33) was absent (C -1 = 0). 
Matching of the fields yields 

rI (t/o 1 - C/3) [ch(« n (i2 w - #)) + re n i2sh(K n (i2 w - R)j 



a sh(K 1 R)sh(K 1I (R w - i?)) + ch(/t I J R)ch(/t II ( J R w - R)) ' 



(34) 



- [ciuyg) - sh^y^)] 



a sh(/t I J R)sh(/t II ( J R w - #)) + ch(/t I J R)ch(fi; II ( J Rw - #)) 
where we introduced notation a = kP/k 1 . 

The charge in the sphere of current radius r < R is given by 

Q(r) = rch(A) ^1 - < 0, (36) 

being, thereby, negative, since [/" > U\ and V^ 1 < 0. This negative charge is completely 
screened by positive charge induced in the region R < r < i?w 

Then we calculate contribution to the thermodynamic potential (effective energy) of the 
Wigner-Seitz cell per droplet volume. We start with the proper electric field energy term; 



~ _ ~i , ~ii 
ey — e v + e v , 



3 f R {VV l ) 2 A 2j 



AttR 3 Jo 8ne 2 

(1 + aoeth(a 1 Q) 2 (-|th 2 e + | - |th 2 e + f th£) 

~ e (a the • th(aiO + 1) 2 ' ( } 



being expressed in dimensionless units 



Qi- ^5 , / ft- ^ , (38) 

where we used eqs. (18), (28), (34). With the help of eqs. (31), (35), from eq. (18) we find 



y 4vri? 3 ii? 87re 2 

, 2 



(l - |thg) (l - |a aie 2 (l - thVfl) + |aogth(aiO) 

" e 2 (MK • th( ai + 1) 2 ' (39) 
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In order to explicitly calculate correlation terms we introduce the quantity 5V l = V l {r) + 
li e for r < R, thus taking V} = —/i e . Averaging (22) over the droplet volume, with the help 
of (28), (34), we obtain 

_ 3 r\ , ^ p (i + ^th^O) 2 U t M 



2ft (l + aoftMaxQ) _ ft 

+ e (aothe • th( ai + 1) 1 4 U 3' (4Uj 

In the hadron phase introducing SV U = V ll {r) — V ll (r = i? w ), where we used V} 1 = V ll (r = 
R w ) ~ —ii e , with the help of eqs. (22), (31), (33), (35), we obtain 

^ = T^L 4*r 2 dru&* (41) 



, 2 



fta (l - |the) [|th(aiO + 3ai(l - thVQ) - |th(a 1 0/ch(a 1 Q 



6 (a th£ • th(aiO + l) 2 

One can see that u)" r — > 0, if a.\ — > 0, and also in the case a — 0. We could also use other 
values for V} and K r n , e.g. we have checked that using general eq. (21) with V} = V l {ti) and 
V^ 1 = V Il (Rw) leads to the very same result. 

In our dimensionless units the total quark plus hadron surface contribution to the energy 
per droplet volume renders 

es/A) = ft/£, Pi = 3i£<t/Po, (42) 

see (38), and we used that es = 3a/R. Coefficients ft, ft are evaluated with the help of 
eqs. (29), (38) and (42). For the above used quantities fj, e ~ 170 MeV, /i n ~ 1020 MeV, 
a c ~ 0.4 and m s ~ 120 -h 150 MeV we estimate ft ~ 1.6m*. Thus, with the value a ~ 1.3m^ 
we obtain ft ~ 0.7, whereas with a ~ 10 MeV/fm 2 ~ 0.14m^ we would get ft ~ 0.08. 

Coulomb solution for the case of a tiny quark fraction volume is obtained, if we first put 
ai — > oo, and then expand the terms e v + Zy + es in £ C 1. Thus, we recover the Coulomb 
plus surface energy per droplet volume 



(43) 
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where partial contributions correspond to the terms e v , €y and es- Both the correlation 
terms u)\ or oc £ 4 and u;" r oc £ 3 can be dropped in the Coulomb limit, for droplets of a tiny- 
size £ <C 1. 

Function ec,s has the minimum at £ = £ m = (15/31/4) 1 / 3 , corresponding to the optimal 
size of the unscreened droplet. Coulomb solution is reproduced only for ^ m < 1, whereas 
with above estimate fli £ 0.1 we always get £ m ~ 1. On the other hand, for (5\ <C 0.1 we 
would obtain k 1 > m n , corresponding to unrealistically small droplet size R < l/m n . Thus, 
we conclude that pure Coulomb solution is never realized within mixed phase. 

In the limit cki£ ^> 1, £ ^> 1, corresponding to the single large size drop, from (37), (39), 
(40) and (41) we find that all the terms contribute to the surface energy density (oc £ _1 ) 
and, therefore, the electric field effects can be treated with the help of an effective surface 
tension. The full surface tension o"to t her then renders 



spher „,„ „ \I Po a o[ a O + 4/3] 

0tot =a + a v = a- X D — - . (44) 

3(1 + a ) z 



The first a term is the contribution of the strong interaction, and the second negative 
term is the contribution of the electric field effects, which depends largely on the values of 
parameters. For /x e = 170 MeV, fj, n = 1020 MeV, m s = 150 MeV, a c = 0.4, we estimate 
the contribution to the surface tension from the electric effects as ay ~ —70 MeV/fm 2 . 

In Figure we demonstrate dependence of the contribution of inhomogeneous charge 
distributions to the total thermodynamic potential per droplet volume (solid curves), 
SuJtot/Po = (ey + &lor + <^oor + / Po, given by the sum of partial contributions (37), (39), 
(40), (41) and (42), for the case of a single droplet, as function of the droplet size £ for three 
values of a at fixed value f3\ (each panel). The curves labeled by "C" demonstrate the 
Coulomb solution ec,s/A), determined by eq.(43). Each Coulomb curve has a pronounced 
minimum at £ = £ c oc (3l /3 . For £ > £ c the Coulomb curve shows quadratic growth deviating 
drastically for £ > 1 from the solid curve. For fi 1 < 0.01, minimum points of the Coulomb 
curves £c deviate little from the minima of the solid ones. Only for such small values of 
Pi and £ m we recover the Coulomb limit! However, one may obtain such small values of (3\ 
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only for tiny values of surface tension and very large values of neutron chemical potential. 
With increase of the latter, the Debye parameter k 1 is also increased and the droplet radius 
R = Cm/^ 1 is proved to be essentially smaller than l/m n . For larger values of Pi deviation 
between the minima of total and Coulomb solutions is proved to be pronounced. All the 
solid curves converge to 1/3. Large £ asymptotic of solid curves is oc l/£ being interpreted 
as the surface energy term, characterizing by a significantly smaller value of surface tension 
(44) than that determined only by strong interaction. We see that for a > | ay | that cor- 
responds to Pi > Pic ( according to our Figure Pi c ~ 0.5 ) the structured mixed phase is 
proved to be prohibited, since necessary condition of its existence (presence of minimum 
in the droplet size) is not satisfied, whereas with the Coulomb solution necessary condition 
is always performed. Large size drops are realized within the mixed phase if a <| oy |, 
| a + oy |<C| oy | and within Maxwell construction if a >| ay |. The "e.m." curve shown in 
each panel demonstrates the contribution e em /Po = (ey + £s)/Po, ignoring correlation terms 
for a — 1. This quantity e em is the counterpart of the Coulomb solution ec,s when the 
charge screening effect is taken into account. We see that the minimum at the "e.m." curve 
disappears already at Pi > 0.03. Difference between the solid and "e.m" curves shows the 
important contribution of correlation energy in the hadron-quark structured mixed phase. 

Dependence of the curves on the ratio of the screening lengths «o is also rather pro- 
nounced, whereas it was completely absent for the Coulomb solution. Our calculations also 
show that dependences of 5u tot /Po on the volume fraction / are very weak in the whole 
range of available values. 

Summarizing, in discussion of possibility of presence of structured mixed phase at first 
order phase transitions in multi-component systems of charged particles we consistently 
incorporated effects of the charge screening. As an example, our formalism was applied 
to the hadron-quark structured mixed phase. We demonstrated that the charge screening 
effect should greatly modify the description of the mixed phase, changing its parameters and 
affecting the possibility of existence. In absence of the mixed phase our charged distributions 
describe the boundary layer between two separated phases existing within the double-tangent 
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(Maxwell) construction. Consideration of non-spherical droplets (rods and slabs) does not 
change our conclusions. Further discussions are given in another paper [8]. 
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Figure: Contribution to the effective energy due to inhomogeneous charge distributions 
per droplet volume in case of a single droplet versus scaled droplet radius. 
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